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Abstract 

A low-energy 77 collider has been discussed in the context of a testbed for a 77 
interaction region at the Next Linear Collider (NLC). We consider the production of 
heavy mesons at such a testbed using Compton-backscattered photons and demon- 
strate that their production rivals or exceeds those by BELLE, BABAR or LEP 
where they are produced indirectly via virtual 77 luminosities. 
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1 Introduction 



Quantum Chromodynamics (QCD) is a gauged theory of the strong interac- 
tions. Exquisitely predictive in the perturbative high-energy regime, at present 
a description of the physical states is beyond the reach of the theory. Me- 
son spectroscopy allows a more complete understanding of QCD in this non- 
perturbative realm. 



Heavy meson spectroscopy provides unique insight into the underlying the- 
ory without the additional complications associated with the the lighter (and 
more relativistic) mesons. Heavy mesons have been studied at e + e~ colliders 
for nearly 30 years, however, only selective states compatible with the quantum 
numbers of the e + e~ beams are produced directly: other states are populated 
via subsequent hadronic or electromagnetic decays. Often, as is the case with 
the 7/& (nS) mesons, the relevant branching ratios are very small. Photon in- 
teractions provides a mechanism by which some of these inaccessible mesons 
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are produced directly. In a way that is complementary to e + e~ interactions, 
77 interactions produce particles having final states with C of +1, where C is 
the charge conjugation number. 

Participating photons can be either virtual or real. Virtual photons are an 
irreducible description of charged particle beams: all e + e~ colliders generate 
virtual 77 luminosity - in the case of PEP-II or KeK-II the virtual luminos- 
ity can rival the e + e~ luminosity. An alternate mechanism for producing 77 
luminosity is with real photons. Typically, low energy photons are made to 
Compton-backscatter off of an energetic charged particle beam to produce 
highly boosted photons. Opposing photon beams are kinematically focused 
and interact to produce real 77 luminosity. Although this production mech- 
anism requires the added complication of introducing a laser and associated 
optics into the interaction region, the substantial benefits include control over 
the 77 luminosity profile, polarization and magnitude. Technology advances 
in optics and lasers have progressed to the point that one could construct a 
low energy 77 testbed as a proof-of-principle for the NLC, where it would 
be possible to singly produce the Higgs boson (whereas it must be produced 
in pairs at e + e~ colliders). The potential this offers for doing precision Higgs 
studies is a major driving force behind the impetus for a 77 option at the NLC. 
For the purposes of this study, however, we restrict our attention to photons 
with a center of mass energy of 60 GeV, which is sufficient for an engineer- 
ing demonstration. One possible site for the testbed is the Stanford Linear 
Collider, where polarized electron beams meet opposing positrons beams at 
center-of-mass energies up to the mass of the Z . In this proceeding we demon- 
strate that with suitable choice of laser and optics parameters the resulting 77 
luminosity can be made to exceed the virtual 77 luminosity from e + e~ collid- 
ers. This luminosity can thus serve as the basis for a physics program at the 
testbed facility, e.g., heavy quarkonia studies. In subsequent sections we com- 
pare and contrast real and virtual 77 luminosities and discuss the implications 
for heavy quark meson spectroscopy. 



2 Real photon luminosity 



Depending on the energy transfer scale the structure of an otherwise feature- 
less photon can be resolved, resulting in an irreducible hadron and leptonic 
background that would complicate precision Higgs measurements. Hence, a 
primary goal of a 77 testbed would be to measure the 77 luminosity and com- 
pare the result with simulation. In this study Compton-backscattered lumi- 
nosity was simulated with CAIN 2.1e (0) using the laser and beam parameters 
listed in Table 1, where the electron and laser polarizations refer to circular 
polarization. Beam parameters in this table reflect NLC-like expectations 
about achievable beam sizes, while the laser parameters were chosen to maxi- 
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Beam 


Laser 


parameter 


value 


parameter 


value 


beam energy 


30 GeV 


frequency 


120 Hz 


charge 


4xl0 lu 


pulse energy 


2.0 J 




1.47 jim. 


A 


1.05 /J,m 


Gy 


0.0552 jum 


laser polarization 


± 1 




0.1 mm 






X 


8.0 mm 






Py 


0.1 mm 






e~ polarization 


0.8 






e + polarization 










Table 1 

Beam and laser parameters used to characterize the real photon luminosity at LINX. 



mize the Compton-backscattering rate. (Throughout this proceeding we refer 
to a hypothetical accelerator whose beam parameters match those of Table 
1 as LINX.) Figure 1 shows the resulting simulated Compton-backscattered 
electron and resultant photon spectra at the LINX conversion point. With 
nearly 2xl0 10 photons per electron the conversion is very efficient: more than 
98% of the electrons undergo Compton-scattering. The preponderance of pho- 
tons also causes significant degradation of the single-interaction photon and 
electron Compton edges around 10 and 20 GeV, respectively. This degradation 
is due to multiple Compton interactions as the electron beam traverses the 
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Fig. 1. A snapshot of the right-going photons and electrons at the conversion point. 
The ratio of backscattered photons to initial electrons is approximately 4:1. Note 
the different axis scalings for photon and electrons. 
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laser pulse. This picture is quite distinct from the sharp Compton edges seen 
in simulations (not shown) where the energy of the laser pulse is considerably 
lower and single interactions dominate. 

77 luminosity is produced when the left- and right-going photons traverse the 
several millimeters from their respective conversion points to the interaction 
point. Figure 2 shows the resulting simulated Compton-backscattered 77 lu- 
minosity at the interaction point from CAIN. For this analysis the luminosity 
is divided into 50 equally spaced energy bins spanning the range < y/s < 
62.4 GeV. Multiple Compton interactions give rise to significant luminosity 




5 10 15 20 25 30 

E (GeV) 



Fig. 2. 77 luminosity for the beam and laser parameters in Table 1. The total 
integrated (instantaneous) luminosity is 2.971 x 10 32 cm" 2 s _1 

above 20 GeV, whereas for lower laser energies (where single interactions dom- 
inate) the luminosity is essentially zero for energies above twice the Compton 
edge. The area under the curve in Fig. 2, or the instantaneous 77 luminosity, 
is 2.971 xl0 32 cm~ 2 s _1 . By way of comparison, the maximum instantaneous 
e + e~ luminosity at LEP never exceeded 1.0xl0 32 cm~ 2 s _1 over its 10+ years 
of operation. 



3 Virtual photon luminosity 

To assess the relative merits of a real photon collider we have considered 
three colliders whose large e + e~ luminosities would suggest a large virtual 
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77 luminosity: PEP-II (2), KEK-B (0) and LEP (Q). PEP-II and KEK-B are 
asymmetric e + e~ colliders running at a center of mass energy of 10.58 GeV, 
while LEP was a symmetric e + e~ machine that scanned over the energy range 
91.2 <^/s < 206 GeV until decommissioned in late 2000 to make way for the 
Large Hadron Collider (LHC). 

The classical treatment of virtual photons exploits the equivalence between the 
fields of energetic charged particles and pulses of electromagnetic radiation. 
In the equivalent photon approximation (EPA) a pulse of electromagnetic 
radiation is transformed into a frequency spectrum of virtual photons. The 
interaction of virtual photons from colliding charged particle beams gives rise 
to virtual photon luminosity. The Weizsacker Williams approximation to the 
EPA restricts itself to transverse modes of the virtual photon while integrating 
over the (unobserved) momentum transfer between the electron and the virtual 
photon to yield (0) 



dN a 



dz 27T z 



\-z+^)\o g { q ^)-2ml(4 

, / Qmin \ Qmin Qmax I 



(1) 



where z is defined as E y /E beam , m e the electron (or positron) mass, q 2 min and 
Imax are the minimum and maximum photon energies, respectively, and dN/ dz 
is the number of photons in an energy interval dz. Both g^ in and q^ax depend 
on z in the following manner 



2 m l z2 (2) 

Hmm \ — z ^ ' 

(Lax = (1 - Z ) E beam 2 i 1 ~ ™S 9 max ) (3) 



The dependence of q^ax on @max in Eq. 3 merits comment. As 9 max increases 
the opportunity for detecting electrons diminishes due to limited detector 
acceptance, but it is precisely these unobserved electrons that contribute to Eq. 
1. The appearance of 6 max the above expression requires that specific detector 
configurations be taken into account in the WW analyses. BABAR (6) and 
BELLE (0) detectors surround the interaction regions at PEP-II and KEK-B, 
respectively. We take the OPAL (0) detector as the prototypical LEP detector, 
while the acceptances of the LINX detector were assumed to be those of the 





BABAR 


BELLE 


OPAL 


LINX 


Energy (GeV) 


3.1/9.0 


3.5/8.0 


91-206/91-206 


30/30 


9 max (mrad) 


523.6/295.6 


715.6/401.4 


33/33 


36/36 



Table 2 

e + /e~ beam energies and forward/backward acceptances for the BABAR, BELLE, 
OPAL and LINX detectors. 
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year 


BABAR 


BELLE 


OPAL 


peak 1 


int 2 


peak 


int 


peak 


int 


1990 


- 


- 


- 


- 


0.0110 3 


0.0121 


1991 


- 


- 


- 


- 


0.0110 3 


0.0189 


1992 


- 


- 


- 


- 


0.0110 3 


0.0286 


1993 


- 


- 


- 


- 


0.0190 


0.0400 


1994 


_ 


_ 




_ 


0.0231 


0.0645 


1995 


- 


- 


- 


- 


0.0341 


0.0461 


1996 


- 


- 


- 


- 


0.0356 


0.0247 


1997 


- 


- 


- 


- 


0.0570 


0.0734 


1998 


- 


- 


- 


- 


0.0999 


0.1997 


1999 


1 


1.62 


2.05 4 


0.287 


0.1069 


0.2537 


2000 


2 


23.76 


2.05 


10.940 


0.0675 


0.2331 


2001 


4.4 


40.05 


5.171 


36.257 






2002 


5 


31.32 


8.256 


54.181 






2003 


7.5 


56.71 


11.305 


77.232 






2004 


7.5 5 


4.63 


11.305 5 


6.827 






Total 




158.09 




185.724 




0.9948 



Table 3 

Instantaneous and integrated e + e~ luminosities as a function of year. For BABAR 
and BELLE the 2004 luminosities are through 01/04. 

Stanford Linear Detector (0). Table 2 lists the beam energies and forward and 
backward acceptances seen by the four detectors under consideration. The 
forward and backward acceptances of BABAR and BELLE are unequal due 
to the asymmetric beam energies. 

In order to compare real and virtual 77 luminosities yearly e + e~ luminosity 
records for BABAR, BELLE and OPAL were compiled. Table 3 summarizes 
the peak and integrated e + e~ luminosities seen by each detector since initia- 
tion of their respective operations. 



4 Results 



To simulate Weizsacker Williams luminosity the physics event generator PAN- 
DORA 2.1 (p) was used. Since Eq. 1 increases without limit as z— >0, it was 
necessary to introduce a low-energy cutoff to regularize the result. An energy 
cutoff of 1 MeV was employed in all subsequent analyses. The choice of cut- 
off obviously affects the integrated virtual luminosity, but has no discernible 

1 In units of nb _1 s _1 

2 In units of fb _1 

3 LEP peak luminosities assumed for 1990-1992. 

4 BELLE peak luminosity not available for 1999. 

5 BELLE and BABAR peak luminosities for 2004 taken from 2003. 
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affect on the virtual photon spectra above the energy range of interest (~ 3 
GeV), nor on the luminosity integrated cross sections discussed below. Virtual 
77 luminosity was calculated for the PEP-II, KEK-B and LEP beams using 
Pandora in two steps. First, a random number corresponding to the energy 
of a virtual in one beam was generated. Associated with this photon energy 
is a specific value of dN/dz from the Weizsacker Williams approximation. 
This same procedure was repeated for the opposing beam. The product of the 
dN/dz from each opposing beam gives the relative virtual 77 luminosity per 
unit e + e~ luminosity, which for PEP-II was determined to be 0.19. Instanta- 
neous virtual 77 luminosities were derived by multiplying the relative virtual 
77 luminosities by the instantaneous e + e~ luminosities from Table 3. Although 
the table contains yearly peak luminosities (when available), for our purpose 
the peak luminosity with the greatest magnitude was used. A comparison with 
the LINX instantaneous 77 luminosity is shown in the lower graph in Fig. 3. 
In both graphs in this figure and in the analyses described below the BELLE 
and BABAR luminosities have been summed together, while the LEP lumi- 
nosity is that seen by the OPAL detector, but multiplied by four to account 
for the four LEP detectors. The characteristic rise inherent in the Weizsacker 




Fig. 3. The instantaneous (lower) and integrated (upper) virtual Weiz- 
sacker- Williams luminosities seen by BELLE+BABAR and LEP compared to the 
real luminosity seen at LINX. All luminosities are in 1.3 GeV bins. 

Williams description of virtual photons is seen in the BELLE+BABAR and 
LEP plots, whereas the rise in the LINX luminosity at low energies stems in- 
stead from bremsstrahlung and multiple scattering. Using the relative 77 lumi- 
nosities above and the peak luminosity data from Table 3, one arrives at peak 
instantaneous luminosities of 3.479 x 10 33 cm~ 2 s _1 and 3.121 Xl0 32 cm~ 2 s _1 for 
BABAR+BELLE and LEP, respectively. By comparison, the LINX instanta- 
neous luminosity is 2.971 xl0 32 cm -2 s _1 . 
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BABAR+BELLE 


LEP 


LINX 


norm 


0.19 


0.52-0.70 


1.0 


inst (10 32 cm~ 2 s -1 ) 


34.79 


3.121 


2.971 


intr (10 39 cm _2 s) 


63.86 


2.60 


2.971 



Table 4 

Normalized, instantaneous and integrated luminosities from BABAR+BELLE, LEP 
and LINX. One Snowmass year of operation (10 7 s) was assumed for the LINX 
integrated luminosity. 

A similar approach was used in producing the integrated luminosities in the 
top graph in Fig. 3, but here the relative virtual 77 luminosities are multiplied 
by the yearly integrated luminosities and summed over relevant years to form 
a total integrated luminosity. For instance, in the case of LEP there are 11 
contributions to the sum (years 1990 to 2000) - each year with its own relative 
and integrated luminosity to account for the changing beam energy with time. 
A summary of the results for the normalized, instantaneous and integrated 
luminosities is contained in Table 4. Taken together, Table 4 and the top 
graph in Fig. 3 is one of the our main findings: although the integrated virtual 
77 luminosity from BELLE+BABAR greatly exceeds the integrated real 77 
luminosity expected from LINX, due to the sharply decreasing nature of the 
Weizsacker- Williams spectrum this dominance disappears beyond 3.0 GeV. 
As will be seen in the next section, since much of the physics of interest lies 
above this energy, real 77 luminosity offers a very viable approach to heavy 
meson spectroscopy. 



5 Physics Opportunities 



One goal of heavy quark spectroscopy is to test Quantum Chromodynamics 
predictions of the masses, widths and quantum numbers of these mesons by 
comparing them with their experimental values. Production of pseudoscalar 
mesons is a particularly attractive physics goal for a low-energy 77 collider 
since the mesons quantum numbers J = + ~ preclude them from being 
produced directly at an e + e~ collider. Although the lowest radial excitations 
of the charmonium system has been thoroughly mapped out (with the excep- 
tion of the singlet h(lp) state), spectroscopy of the second radial excitation 
has numerous holes. Even where there have been recent measurements, as in 
the case of the pseudoscalar i] c (2S), disagreement persists as to the excitation 
(mass). Conflict between experiment and theory on the order of 50 MeV in 
mass of the 1] C (2S) is the source of some consternation in the theoretical com- 
munity (JlOj). Enhanced statistics would further allow for the determination of 
unknown branching ratios. The situation is more severe in the bottom system: 
none of the pseudoscalar bottom mesons have been definitively observed ((lit). 
Bottom pseudoscalar mesons apparently are more difficult to produce at e + e~ 
colliders: their heavier heavier masses imply a more non-relativistic wave func- 
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meson 


mass 


Ftot 


V 

i 77 


%(is) 


2.979 


0.0161 


7.4xl0~ 6 


r? c (2S) 


3.654 


0.0161 


7.4xl0~ 6 


XcO 


3.415 


0.0107 


2.6xl0~ 6 


%(1S) 


9.3 


0.014 


5.3xl0~ 7 



Table 5 

Masses and widths of the mesons in GeV. 

tion. This, in turn, enhances the selection rule which in the non-relativistic 
limit forbids E{1) transitions between the T(nS)^^(nS). 

Using the integrated luminosities described in the previous section we have 
studied the production of pseudoscalar and other suitable heavy charm and 
bottom mesons at a low-energy 77 testbed and have compared these results to 
what can be achieved through virtual 77 luminosity. Meson simulations were 
again carried out using the PANDORA event generator. Meson production 
in PANDORA is modeled using the Breit-Wigner approach; necessary modi- 
fications were made to include the widths and masses of the various mesons 
that are summarized in Table 5. To model meson production at LINX CAIN 
2.1 files were imported into PANDORA. These CAIN luminosity files have 
10000 entries representing the partitioning of the energy range < y/s < 62 A 
GeV into 50 equally spaced energy bins of 1.248 GeV each for E 7l and E 72 . 
There are four sets of 2500 entries corresponding to the four helicity combina- 
tions possible when two photons interact. Each of the 10000 entries contains 
E 7l and E 72 , luminosity and overall helicity. PANDORA normalizes the lu- 
minosity so that the sum over all 10000 entries is 1.0. Meson production at 
BABAR+BELLE and LEP was modeled using the luminosity class described 
in Section 3. 



PANDORA can be made to output luminosity integrated cross sections for 
the process of interest, defined as 



/ 



r)T 

m *dE (4) 



In practice, the integral is replaced by a finite sum sampled from a uniformly 
distributed E 7 i and E 72 plane. The utility of the luminosity integrated cross 
section is that when multiplied by the total integrated luminosity the re- 
sult is the number of events for the process under consideration. Table 
6 lists the luminosity integrated cross sections and total number of events 
produced for each meson in the previous table. Note that unlike the con- 
ventional definition of cross section, the luminosity integrated cross section 
depends on details of the environment in which the meson is produced: it 
is lowest for BABAR+BELLE, for which the integrated luminosity is great- 
est. From this same table it is seen that number of r] c (lS) events produced 
at BABAR+BELLE and LINX are approximately equal. Figure 3 provides 
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a consistency check on this result. In the upper plot in that figure the in- 
tegrated luminosities of BABAR+BELLE and that of LINX are seen to be 
approximately equal at an energy corresponding to the mass of the T] C (1S). 
Beyond this crossover point, a greater number of heavier mesons are produced 
at the LINX facility than at BELLE+BABAR or LEP. For the purposes of 
bottom spectroscopy the LINX scenario is clearly superior, with the poten- 
tial for producing an order of magnitude or more //&(1S) mesons than either 
BELLE+BABAR or LEP. 

In order to simulate experimental reconstruction of these mesons at a 77 
testbed it is necessary to specify the decay mode, as well as potential back- 
grounds. Decay of heavy mesons into pp pairs was chosen due to the simplicity 
of the event reconstruction, the relatively large branching ratios and knowledge 
about the competing background process. The dominant background 77 — >pp 
arises when both photons are resolved into hadrons. This process has been 
well-studied for insight into the photon structure function and resolution of 
the ongoing debate concerning diquark versus independent quark models of 
the hadron. Theoretical and experimental (flih efforts have gone into de- 
termining the cross section for this process. In particular, the latest BELLE 
results describe a differential cross section (in nb/sr) for this process of the 
form 

da - 1 i 5 x 10 V 15 ( 1 + C ° s2 ^ (5) 
dEd(cos6) " 1 - 15xlUff h-cos^J (5j 



for I cos(#)| < 0.6 and ^/s > 2.75 GeV. This result is consistent with the di- 
quark model, in contrast to the independent quark model predicts an isotropic 
angular distribution. The normalization factor that appears in Eq. 5 was de- 
rived directly from BELLE data. A separate Pandora class was written to 
generate the 77 — »pp background; when convolved with the LINX luminosity 
described above a luminosity integrated cross section subject to the above 
kinematic limits of 0.0058 nb is obtained, which corresponds to 17337 events 
over a Snowmass year of operation. Similarly, to compute the number recon- 
structed meson events it was necessary to multiply the raw number of meson 
events from the last row of Table 6 by the branching ratio for decay into 
pp pairs. Table 7 summarizes the number of events expected from this proce- 



meson 


BABAR+BELLE 


LEP 


LINX 


a 


events 


a 


events 


a 


events 


Vc(lS) 


3.19e+4 


1.08e+7 


1.72e+6 


1.98e+6 


3.34e+6 


9.91e+6 


r? c (2S) 


1.25e+4 


4.25e+6 


8.43e+5 


9.74e+5 


1.97e+6 


5.84e+6 


XcO 


5.23e+3 


1.78e+6 


3.32e+5 


3.83e+5 


7.66e+5 


2.28e+6 


%(1S) 


3.05e-l 


1.09e+2 


3.25e+2 


3.89e+2 


1.07e+3 


3.17e+3 



Table 6 

The luminosity integrated cross sections and events produced for the three machine 
scenarios considered. 
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dure. The Phythia interface to Pandora allows events to be written in StdHeP 



meson 


BR(pp) 


PP 


Reconstructed pp events 


%(1S) 


0.0012 


1.19e+4 


7.09e+3 


Vc(2S) 


0.00049 


2.80e+3 


1.69e+3 


XcO 


0.00024 


5.47e+2 


3.32e+2 


%(1S) 


0.0012 


4.00e+0 


3.00e+0 


77 ^PP 




1.73e+4 


1.73e+4 


Total 




3.26e+4 


2.64e+4 



Table 7 

Branching ratios into pp final states and total number of signal and background 
events expected for one Snowmass year of running at a 77 testbed. 



format (14), which in turn are read in by the LCDROOT ffl5h - a detector 
simulation and analysis framework based on ROOT |l6h. In this analysis the 
pp final states and background were propagated through a small version of a 
hypothetical NLC detector whose electromagnetic and hadronic calorimeters 
reside in 5 T magnetic field. Details of the detector (fl7f ) are not crucial to 
the results. Passage of pp through the detector simulation results in momenta 
smearing. In the hadronic calorimeter the momenta is inferred. For this anal- 
ysis only inclusive decays to pp are considered, thus all particles reaching the 
hadronic calorimeters were assumed to have the mass of a proton. With this 
information and the total energy measured in the electromagnetic calorimeter 
the four vector is constructed and the invariant mass is calculated. Figure 4 is 




Fig. 4. Reconstructed meson and background events expected from a Snowmass year 
of operation at LINX. Lower and upper plots are identical except for the scaling of 
the axes. Cuts of | cos(#)| < 0.6 and yfs > 2.75 GeV were applied to all data. 

a reconstruction of the invariant mass of events leaving a track in the hadronic 
calorimeter subject to the restriction that | cos(#)| < 0.6 and yfs > 2.75 GeV. 
The total number of reconstructed events (26,494) found in the last column of 
7 is significantly fewer than the 32,527 StdHeP events input into LCDROOT. 
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This is an artifact stemming from the fact that kinematic cuts were not im- 
posed on the meson events within PANDORA but rather inside LCDROOT. 
For the isotropically decaying mesons the fraction lost is consistent with the 
restricted phase space imposed by the angular cuts. The meson events clearly 
stand out above the pp background, although additional analysis is needed 
to determine whether it is possible to separate out Xco and r? c (2S) events. A 
thorough description of these results, including a complete treatment of the 
non-isotropic decay of the Xc2, will appear elsewhere. 



6 Conclusions 

A 77 interaction region has been proposed for the NLC using Compton- 
backscattered photons. Technical challenges integrating the optical, laser and 
mechanical subsystems suggest the desirability of a low-energy testbed. In 
this paper we address one aspect of a physics program that could accompany 
such a test bed - heavy meson production - and compare it with similar me- 
son production expected where virtual photons are the means of production. 
We find that after a single year of operation the number of charmed mesons 
seen at a LINX facility equals or exceeds the number generated at either 
BABAR+BELLE or LEP to-date. LINX is especially dominant when pro- 
ducing bottom mesons, offering the exciting possibility of discovering ??&(nS) 
mesons before the LHC begins delivering data toward the end of this decade. 
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